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laminates  (0)48  6  mm  in  thickness. '  The  sprcimen  was  of 
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2.2  Fatigue  Test 

The  test  was  carried  out  with  an  electric- 
hydraulic-servo-type  testing  machine  with  capacities  of 
80  N  and  9.8  KN  connected  to  a  computor  under  Ae 
conditions  of  a  constant  stress  ratio  (ratio  of  the  min¬ 
imum  load  to  the  maximum)  and  a  constant  rate  of 
change  in  the  stress  intensity  factor  (l/K)dK/da.^  The 
environment  for  testing  chosen  was  in  air  at  23°C  and  50 
percent  RH  and  in  water  at  23°C  and  50°C;  the  rate  of 
repetition  used  was  lOHz  for  the  standard  and  2  Hz 
partly  for  comparison,  in  load-reduction  tests  and  load¬ 
raising  tests,  respectively.  The  fatigue  test  in  water  was 
carried  out  with  a  trial-manufactured  test  machine  incor¬ 
porating  a  horizontal  actuator,  with  the  end  point  of  the 
crack  alone  dipped  in  water.  The  length  of  the  crack  was 
estimated  from  the  compliance  characteristics  in  the 
course  of  tests.  Here,  the  relationship  between  the  stress 
intensity  factor  K  and  the  energy  release  rate  G  was 
computed  using  the  following  formula^ 

K  (MPam*'^)  -  0.144  x  square  root  of  G{N/m) 

2.3  Water  Absorption  Test,  Provisionary  Processing  of 
Test  Specimens  in  Water  Environment 


[Article  by  Masaki  Hojo  of  the  Product  Science  Labora¬ 
tories;  Kiyoshi  Kenmochi  of  the  same  laboratories;  Kei- 
suke  Tanaka  of  the  Engineering  Department,  Kyoto 
University;  Kazutoshi  Endo  of  Lion  Co.,  Ltd.:  “Effect  of 
Water  Environment  on  Propagation  of  Delamination 
Fatigue  Cracks  in  CF/PEEK  Laminates”] 

[Text]  CFRP  laminates  have  been  applied  in  primary 
stage  structural  material  recently  and  an  assessment  of 
its  damage  tolerance  in  such  actual  conditions  has  grown 
in  importance.  The  authors  have  reported  in  this 
connection'  that  a  CFRP  involving  a  thermoplastic  resin 
PEEK  as  the  matrix,  which  has  been  developed  recently, 
displays  a  notably  improved  toughness  and  its  laminates 
have  a  resistance  against  delamination  fatigue  crack 
propagation  which  is  more  than  double  conventional 
brittle  CFRP  laminates  of  epoxy  derivatives.  Neverthe¬ 
less,  the  resistance  of  the  material  against  the  environ¬ 
mental  effects,  a  property  on  which  much  hope  has  been 
placed,  has  so  far  been  scarcely  investigated. 

This  article  deals  with  the  mode  I  delamination  fatigue 
crack  propagation  behavior  to  investigate  the  effects  of 
water  environments  on  CF/PEEK  laminates  and  time 
dependence  of  the  material,  and  their  relevant  mecha¬ 
nisms. 

2.  Experiment  Method 

The  material  used  was  the  APC-2  prepreg  laminates  of 
ICI  Inc.,  which  the  authors  had  ICI  mold  into  one-way 


Small  pieces  excised  from  the  specimen  were  immersed 
in  water  at  23'C  and  90'C  to  investigate  their  water 
absorption.  Behavior  in  water  absorption  at  90°C  is 
shown  in  Figure  1.  The  experimental  results  revealed 
that  the  water  saturation  rate  at  0.17  percent  was,  by  a 
factor  of  ten,  lower  than  for  the  CFRPs  of  epoxy  resin, 
but  the  rate  of  diffusion  is  of  the  same  order,  being  very 
slow.’  At  room  temperature,  it  takes  as  long  as  about  30 
years  to  bring  the  specimen  to  the  saturation  or  equilib¬ 
rium  state.  In  the  present  experiment,  therefore,  the  test 
specimen  was  prepared  under  the  two  types  of  condi¬ 
tions  identical  with  those  used  for  the  epoxy-type  CFRP 
to  investigate  the  effects  of  water. 


APC-2  in  water 


20  40  _  60  80  100 

Time.jr  ((hours  ) 

Figure  1.  Relation  Between  Moisture  Content  and  Time 
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(1)  Preliminary  Immersion 

The  specimen  was  immersed  in  water  for  12  hours 
without  loading  before  the  testing.  In  this  case,  water 
diffuses  only  over  the  surface  layer  of  the  specimen,  but 
the  fracture  area  at  the  end  point  of  the  crack  was  so 
small  (30  pm,  low-propagation-rate  region)  that  the 
region  contained  more  than  50  percent  of  the  water 
required  to  maintain  equilibrium.  The  preliminary 
immersion,  therefore,  conceivably  indicates  short-term 
effects  of  water  environments. 

(2)  Water  Control  to  Reach  the  Equilibrium 

The  specimen  was  subjected  to  an  enhanced  water- 
absorption  test  in  water  for  some  200  days  to  reach 
equilibrium.  The  specimen  underwent  no  changes  in  its 
fracture  toughness  (Kjc  of  4.7  MPa  m‘^^  and  G,c  of  1 270 
N)  under  such  a  state.  The  specimen  having  reached  the 
equUibrium  conceivably  represents  a  long-term  effect  of 
water  environments. 


3.  Results  and  Discussions 

3.1  Study  of  the  Governing  Dynamic  Parameters 

The  authors  in  their  previous  report  suggested  an  equiv¬ 
alent  stress  intensity  factor  range,  AK*,,  the  value  of 
which  is  given  by  the  following  formula,  as  a  parameter 
that  governs  crack  propagation  under  different  R  value 
conditions,  and  confirmed  that  the  parameter  was  appli¬ 
cable  also  for  the  case  of  APC-2.  ‘ 

AK.,  -  AK  (1  -  R)-^  =  AK‘-^K„,/ 

where  y  and  1  -  y  denote  the  parameters  indicating  the 
relative  rate  at  which  the  maximal  load  and  the  repeating 
load  contribute  to  the  propagation  [of  the  crack].  In  the 
present  experiment,  as  previously,  y  and  hence  AK^q 
were  computed  from  the  relationship  between  AK  and  (1 
-  R)  at  da/dN  =  10'*°  m/cycle.  The  value  of  y  was 
small— ranging  from  0.17  to  0.36 — which  indicated  that 
the  propagation  of  cracks  was  largely  dictated  by  the 
repeating  loads. 

3.2  Short-Term  Effects  of  Water  Environments 

The  crack  propagation  rate  in  water  at  23°C  and  50”C 
against  AK<,q  for  specimens  having  been  previously  sub¬ 
jected  to  the  preliminary  immersion  in  water  is  plotted 
in  Figure  2.  The  crack-propagation  rate  is  not  dependent 
on  the  stress  ratio,  but  on  the  AK^q  exclusively  in  the 
entire  relevant  range.  The  behavior  of  the  fatigue  crack 
propagation  in  water  at  both  23°C  and  50°C  differed 
from  that  in  air  (dotted  line  in  the  figure)  in  that  there 
was  bending  at  approximately  da/dN  =  5  X  10' as 
represented  in  the  exponential  relation.  The  exponential 
indeces  varied  from  7.1  to  11  and  from  6.2  to  18  for 
water  applied  at  23°C  and  SO'C,  respectively;  in  either 
case,  the  index  was  found  to  be  smaller  in  the  low 
propagation-rate  region  than  in  the  high.  The  lower  limit 


didn’t  seem  to  exist  under  both  conditions,  and  devel¬ 
opment  of  time  dependence  in  the  low  crack- 
propagation-rate  region  was  conceivable.  When  com¬ 
pared  at  the  same  AK^q,  the  crack  propagation  rate  fell  in 
the  region  of  the  higher  propagation  rate;  the  differences 
were  a  fall  to  one-third  at  water  temperature  of  23°C  but 
to  one-half  to  one-tenth  at  50°C.  The  differences  were 
not  notable  in  the  region  of  low  propagation  rate.  In 
contrast,  914°C  laminates  of  the  epoxy-series  had  its 
crack-propagation  rate  lowered  to  one-tenth  in  23''-C 
water  and  its  exponential  index  raised  enormously  in 
50°-C  water.  The  environmental  effects  of  water  on 
APC-2  laminates  were  limited  in  comparison  with  the 
above  epoxy-laminates,  for  which  one  possible  reason  is 
that  the  CFRP  using  an  epoxy  resin  as  the  matrix  absorbs 
water  up  to  2  to  3  percent  while  APC-2  has  a  low  water 
absorption  of  0. 1 7  percent. 


Figure  2.  da/dN  vs.  in  Water  at  23  and  50°C  After 
Prior  Immersion 

3.3  Long-Term  Effects  of  Water  Environment 

Test  specimens  were  immersed  in  23°C  water  where  they 
were  subjected  to  long-term  control  to  reach  the  satura¬ 
tion  limit.  Results  are  shown  in  Figure  3.  For  each  stress 
ratio,  the  curve  in  the  plot  indicating  an  exponential 
relationship  produced  bending  approximately  at  da/dN 
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Fignre  3,  da/dN  vs.  K«,  in  Water  at  23°C  With  Mois¬ 
ture  Conditioning 


-  5  X  10  '°  m/cycle,  as  did  the  curve  for  the  specimen 
treated  with  the  preliminary  immersion  in  water,  with 
the  exponential  index  smaller  in  the  region  of  lower 
crack  propagation  rates. 

The  effects  of  the  stress  ratio,  in  turn,  differed  between 
the  level  above  the  value  of  the  crack  propagation  rate 
da/dN  of  5  X  10  '°  m/cycle  and  the  one  below;  the  crack 
propagation  rate  in  the  low  rate  region  were  dependent 
on  the  AK«,  alone  [and  independent  of  the  stress-ratio 
values]  while  the  propagation  rate  grew  with  an  increase 
of  the  stress  ratio  value  in  the  high-propagation-rate 
region.  For  the  low-propagation-rate  region  at  R  of  both 
of  0.2  and  0.5,  and  for  the  high-propagation-rate  region 
at  R  =  0.2,  the  relationship  between  AK^q  and  da/dN  for 
the  specimen  subjected  to  the  long-term  immersion  in 
water  was  approximately  identical  with  the  one  for  the 
specimen  subjected  to  the  preliminary  immersion  fol¬ 
lowed  by  the  fatigue  test  at  50°C  so  that  the  resistance 
against  the  crack  propagation  for  the  specimen  with  the 
long-term  water-immersion  was  larger  than  the  crack 
propagation  resistance  for  the  specimens  subjected  to 
exposure  to  air  alone  or  to  the  preliminary  immersion 
followed  by  the  fatigue  test  at  23°C.index  were  observed 


In  the  case  of  the  high-crack-propagation-rate  region 
with  R  -  0.5,  however,  the  specimen  with  long-term 
water-immersion  had  a  crack  propagation  rate  about  five 
times  larger  than  for  those  exposed  only  to  air,  indicating 
substantial  deterioration  in  quality.  APC-2  has  great 
durability  in  a  normal  environment  compared  to  CFRP 
of  the  epoxy-series,  which,  when  subjected  to  the  long¬ 
term  controlled  immersion  in  water,  exhibits  deteriora¬ 
tion  in  the  entire  range  regardless  of  the  stress-ratio 
value.  For  example,  the  specimen  with  the  long-term 
immersion  in  water  had  a  low  limit  value,  0.65  times 
than  those  Subjected  to  the  preliminary  immersion  in 
water.^’** 

3.4  Phenomenon  of  Time-Dependence 

The  plot  of  da/dN  vs.  AK*q  representing  an  exponential 
relation  in  Figures  2  and  3  exhibits  a  knicking  of  the 
curve,  suggesting  a  time-dependence  developed  in  the 
low-propagation-rate  region.  A  2-Hz  load-increasing  test, 
therefore,  was  carried  out  in  air  and  in  water  at  23°C  and 
50-0,  with  the  preliminary  immersion.  The  result  is 
shown  in  Figures  4-6  adding  the  results  of  a  relevant 
10-Hz  experiment  in  dotted  lines  for  comparison.  With 
results  for  R  of  0.2  alone  available,  AK^q  was  calculated 


Figure  4.  da/dN  vs.  Keq  in  Air  at  2  Hz 
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using  y  values  at  10  Hz.  No  effects  of  the  rate  of 
repetition  was  notable  in  the  air  as  seen  in  Figure  4, 
whereas  an  increase  in  the  crack  propagation  rate  and  a 
fall  in  the  low-propagation-rate  region  of  the  exponential 
in  water,  with  the  difference  more  pronounced  in  exper¬ 
iments  at  50°C,  as  shown  in  Figures  5  and  6.  In  addition, 
no  knicking  in  the  curve  for  the  exponential  relationship 
was  notable,  the  exponential  index  in  the  entire  crack- 
propagation-rate  region  approaching  that  of  the  low- 
crack-propagation-rate  region  of  the  10  Hz  curve.  A 
comparison  of  the  rates  for  2  Hz  and  10  Hz  in  water  at 
50°C  for  the  range  AK^g  =  1 .5  -  2.0  MPa  m*^^  proved  they 
differ  by  a  factor  of  five,  equivalent  to  the  difference  in 
the  rates  of  repetition.  This  implies  that  the  crack 
propagation  in  this  region  is  dictated  not  by  repetition, 
but  by  time  dependence.  The  knicking  seen  in  the  curves 
of  Figures  2  and  3  may  conceivably  result  from  the 
development  of  time  dependence  in  the  low- 
propagation-rate  region. 

3.5  Observation  of  the  Fracture  Surface 

Figure  7  [not  reproduced]  presents  photographs  taken 
with  a  scanning  electron  microscope  that  show  fracture 
surfaces  produced  in  static  (fracture-toughness)  tests  and 


Figure  5.  da/dN  vs.  Keq  at  2  Hz  in  Water  at  23'C 
After  Prior  Immersion 


in  fatigue  tests  carried  out  either  in  air  or  in  water  at 
23°C,  with  the  specimen  subjected  to  the  preliminary 
immersion  or  the  long-term  controlled  immersion  in 
water  for  the  fatigue  test  in  water  exclusively.  In  (a),  the 
fracture  surface  having  undergone  the  static  test  exhibits 
a  pronounced  plastic  deformation  of  the  resin  and  is 
characterized  by  an  elongation  of  the  resin  in  the  direc¬ 
tion  normal  to  the  surface:  in  (b),  the  fractured  surface 
after  the  fatigue  test  in  air  is  router  indicating  that  the 
damage  extends  deeper  toward  the  inside,  but  no  elon¬ 
gation  of  the  resin  is  notable.  In  specimens  tested  in 
water  at  23°C  with  the  preliminary  immersion  in  water, 
the  fracture  surface  is  similar  to  the  one  in  (b)  for  the 
high-crack-propagation-rate  region,  but  some  elongation 
of  the  resin,  as  shown  in  (c),  is  notable  for  the  low  rate 
region.  For  the  specimens  subjected  to  the  long-term 
controlled  immersion  in  water,  the  elongation  of  the 
resin  is  notable  in  parts  for  even  the  high  rate-region,  and 
most  obvious,  as  shown  in  (d),  in  the  low-rate  region. 
Since  an  elongation  of  resin  more  like  that  for  static 
fracture  surfaces  was  observed  in  the  low- 
crack-propagation-rate  region,  the  mechanism  of  the 
development  of  this  phenomenon  may  conceivably  be 
closely  related  to  the  increased  ductility  of  the  resin. 


Figure  6,  da/dN  vs.  K.,  at  2  Hz  in  Water  at  SO'C 
After  Prior  Immersion 
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Numerous  reports  have  been  given  on  the  enhancement 
by  water  absorption  of  the  creep  deformation  of  polymer 
materials. 
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[Text] 

1.  Introduction 

Three-dimensional  fabric-reinforced  composite  mate¬ 
rials  (3DCM)  may  allow  the  application  of  advanced 
composite  materials  to  extend  to  plastic  materials  of 
great  thickness  and  in  complicated  shapes  since  it  allows 
laminated  composite  materials  (laminated  material)  to 
be  designed  in  a  three  dimensional  form.  However,  the 
three-dimensional  fabric  (3DF),  a  reinforced  3DCM, 
involves  difficulties  in  weaving  and  the  effect  of  three- 
dimensional  arrangement  of  reinforced  fiber  on  the 
properties  of  the  fabric  has  yet  to  be  examined  in  many 
points  with  its  basic  research. 

The  present  report  concerns  a  3DCM,  using  a  3DF  with 
mutually  rectangular  x,  y,  and  z  axes  as  the  reinforcing 
material,  and  examines  its  fracture  toughness  in  compar¬ 
ison  with  the  behavior  of  laminated  materials. 

2.  Material  to  Be  Tested 

Two  types  of  test  materials  were  trial-manufactured;  one 
3DCM  had  filaments  of  24  Kf  with  a  volume  fraction 


(Vf)  for  the  fiber  of  0.40  and  the  other,  those  of  12  Kf 
with  a  Vf  of  0.49,  as  shown  in  Table  1.  Laminated 
materials  with  filaments  1  Kf  and  6  Kf  and  with  Vf  in  the 
range  of  0.34  to  0.66  were  also  trial-manufactured  for 
comparison. 


Table  1.  Raw  Material 


Item 

3DCM 

Laminate 

Fiber 

Type  of  fiber 

Torayca  T300 

Torayca  T300 

Filaments 

6Kfx2  12Kfx2 

lKf6Kf 

Fabric 

Vf  (total) 

0.49  0.40 

0.34  -  0.66 

X  pitch  (mm) 

0.92  2.43  (Thk) 

1.11  mm  2.5 
mm 

Y  pitch  (mm) 

1.05  2.43  (Wld) 

1.11  mm  2.5 
mm 

Z  pitch  (mm) 

3.60  3.00  (Lng) 

Matrix 

Epoxy  resins 

Bisphenol  F  diglycidyl  ethel 

Hardener 

Methyl  hexahydrophthalic  anhy¬ 
dride 

Catalyst 

Tridimethyl  aminothylphenol 

Cure  tempera¬ 
ture 

60“C  (15h)  +  90"C  (2h)  +  120'C 
(2h) 

All  of  these  reinforced  materials  included  epoxy  resin  as 
shown  in  Table  1,  incorporated  by  vacuum  inpregnation 
for  3DCM  and  by  press  molding  for  the  laminated 
materials.  The  conditions  for  hardening  for  both 
molding  methods  were  identical.  The  molded  materials 
of  size  300mmX30mmX  150mm  for  3DCM  and 
llOmmXllOmmXllmm  for  the  laminated  material 
were  cut  to  produce  the  chevron  notch  test  specimens 
(CN  test  specimens)  for  measurement  of  the  work  of 
fracture  (W.O.F.)  as  shown  in  Figure  1.  The  laminated 
material,  as  seen  in  Figure  2,  was  processed  into  test 
tnatfriak  with  two  different  directions  of  lamination; 
one  flatwise  and  one  edgewise.  In  carrying  out  the 
bending  tests,  the  direction  of  the  maximum  density  was 
made  to  match  the  lengthwise  direction.  Table  2  gives 
mechanical  properties  measured  on  these  test  specimens 
to  compare  their  static  and  dynamic  destruction  behav¬ 
iors. 


Figure  1.  Measurement  of  W.O.F.  With  Cheyron-Notch 
Test  Specimen 
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Figure  2.  Shapes  of  the  Test  Specimens  of  a  Flatwise 
and  Edgewise  Laminate 


Table  2.  Comparison  of  Properties  Other  Than  Fracture 
Toughness  Between  3DCM  and  a  Laminated  Board 


Sample 

Com¬ 

pressive 

strength 

(MPa) 

Flexure 

strength 

(MPa) 

Shear 

strength 

(MPa) 

Impact 

energy 

(I) 

3DCM 

350 

454 

162 

11.7 

Laminate  (Flat¬ 
wise) 

- 

323 

25 

5.9 

(Edge-wise) 

- 

333 

80 

2.6 

(Longitudinal) 

228 

- 

- 

- 

3.  Method  of  Assessing  Fracture  Toughness 

W.O.F.  and  the  (strain)  energy  release  rate  (Gr)  were 
computed  from  the  load-displacement  curve  obtained  by 
the  three-point  bending  test  shown  in  Figure  1.  The 
measurements  were  taken  with  a  universal  material¬ 
testing  machine  of  the  Instron  corporation,  Model  1185, 
at  a  crosshead  speed  of  0.5  mm/min.  The  CN  test- 
specimen  was  used  because  it  has  special  advantages 
such  as  its  stable  growth  of  cracks  eliminating  an  intro¬ 
duction  of  initial  cracks  and  that  it  is  suitable  for  thick 
materials,  hereby  providing  the  assessment  of  3DCM 
with  rod-like  material. 

The  W.O.F.  represents  the  amount  of  work  for  a  unit 
cross  section  obtained  on  the  assumption  that  the  total 
eternal  work  applied  to  the  test  specimen  is  consumed 
in  fracture,  and  calculated  by  dividing  the  area  singled 
out  by  the  load-displacement  curve  and  the  axis  of 
abscissa  in  Figure  3  by  an  area  twice  that  of  the  projec¬ 
tion  of  fractured  surface.  Whereas  in  isotropic  material 
the  W.O.F.  is  equivalent  to  Gr,‘  in  anisotropic  material 
Gr  generally  depends  on  the  length  of  the  cracks  pro¬ 
duced,  thus  exhibiting  the  so-called  R  curve,^  and  hence 
some  relevant  information  is  necessary  in  addition  to  the 
average  value  W.O.F.  The  dependence  on  crack  length  of 
Gr,  therefore,  was  computed  from  the  load- 
displacement  curve  by  applying  Bluhm’s  slice-model’ 
and  the  method  of  Munz,  et  ah'*  The  removal  and 
reapplication  of  loads  was  repeated  in  the  figure  with  the 
view  to  obtaining  the  compliance  which  was  necessary 
for  calculating  the  Gr. 


4.  Results  and  Discussion 

Results  of  measurements  for  W.O.F.  are  shown  in  Fig¬ 
ures  4(a)  and  4(b),  which  involves  flatwise  laminated 
materials  and  edgewise  laminated  ones,  respectively,  for 
comparison  with  3DCM,  and  wherein  the  axis  of 
abscissas  represents  the  volume  fraction  of  the  fiber 
arranged  lengthwise  to  the  specimen  bearing  the  load.  It 
can  be  seen  from  these  figures  that  the  flatwise  laminated 
materials  has  higher  W.O.F.  values  than  those  edgewise 
laminated  and  that  the  3DCMs  exhibits  still  higher 
W.O.F.  values.  It  is  also  interesting  to  note  that  the  6Kf 
laminated  material  tends  to  have  larger  values  than  the 


Figure  3.  Plot  of  Stress  vs.  Strain  for  Chevron-Notch 
Specimens 


W.O.F.(KJ/m^)  W.0.F.(KJ/m2) 
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(O) 


(b) 

Figure  4.  W.O.F.  Measurement  by  a  Chevron-Notch 
Three-Point-Bending  Test— (a)  Comparison  between 
3DCM  and  a  flatwise  laminate;  (b)  Comparison  between 
3DCM  and  an  edgewise  laminate 


IKf  ones.  A  notable  difference  of  the  fracture  process  of 
3DCM  from  that  of  the  laminated  materials  is  a  separa¬ 
tion  of  composite  materials  at  interfaces  between  fiber 
bundles  and  the  matrices  accompanied  by  bridgings 
generated  therein.  Figure  5  indicates  development  of  a 
large  number  of  separations,  long  ones  running  normal 
to  the  plane  of  the  fractured  surface,  and  those  running 
parallel  to  the  plane  at  sites  slightly  apart  from  the  plane. 
The  laminated  material  lacks  these  separations  totally 
or,  if  there  are  any  separations,  they  are  very  short. 

Whereas  the  thermal  stress  developed  in  laminated 
materials  in  the  process  of  their  cooling  to  room  temper¬ 
atures  during  their  molding  is  released  in  the  direction  of 
layers,  the  3DCM  has  relevant  strains  built  up  within 
because  of  the  effect  of  the  third  axis  present  and 
resultant  stress  applied  there  may  conceivably  lead  to  the 
long-running  separation.  Development  of  separation  in 
the  interface  even  by  thermal  stress  alone  has  been 
actually  observed  in  some  cases  where  a  slightly  more 
brittle  resin  was  used.  It  is  probable  that  hi^  thermal 
stress  is  generated  internally  even  if  separation  didn’t 
occur. 

Such  long  separation  in  the  bending  test  may  conceiv¬ 
ably  have  resulted  from  sliding  at  interfaces  between  the 
fiber  and  the  matrix  as  proved  by  the  fact  that,  of  the 
three  curves  of  deloading  or  decreasing  loads  in  Figure  3, 
the  one  for  3DCM  alone  does  not  return  to  the  zero-load 
level,  retaining  a  permanent  deformation.  At  the 
advancing  front  of  the  crack,  at  this  juncture,  a  bridging 
of  the  crack  with  fiber  bundles  as  shown  in  Figure  6 
occurs.  As  investigated  and  reported  on  composite  mate¬ 
rials  with  brittle  matrices,  bridging  with  fiber  bundles 
acts  to  block  the  crack-opening  displacement  (COD), 
which  provides  a  highly  probable  mechanism  for  the 
toughness  of  composite  materials.^  Since  the  bridging 
effect  is  more  notable  as  the  interface  separation  grows 
larger,  though  relevant  quantitative  estimation  has  yet  to 
be  carried  out,  the  mechanism  seems  to  explain  in  detail 
the  experimental  results  for  the  toughness  of  3DCM. 


Bridging  of  fiber  bundle  doesn’t  possibly  occur  in  the 
absence  of  cracks;  the  bridging  manifests  its  effect  as  the 
crack  advances.  Figure  7  shows  the  R  curve  effect 


Figure  5.  Damages  Produced  by  W.O.F.-Test  for  3DCM 
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Figure  7.  Plot  of  Strain-Energy-Release-Rate,  G*,  ys. 
Crack  Length 

involved,  with  the  axis  of  ordinate  representing  Gr  and 
the  axis  of  abscissa,  the  length  of  cracks  a,  as  estimated 
by  compliance. 

Though  there  remains  some  room  for  dispute  in  the 
direct  application  of  the  method  of  linear  fracture 


dynamic  for  isotropic  materials,  only  in  the  case  of 
3DCM,  can  cyclic  increase  and  reduction  of  the  Gr  by 
bridging  and  fiber  breakdown  be  observed  clearly.  The 
interval  between  the  Gr  peaks  presents  a  similar  value 
similar  to  the  configuration  of  fiber  bundles,  therefore 
the  energy  created  through  fracture  is  absorbed  in  each 
layer  in  large  quantity. 

The  great  toughness  of  the  3DCM  discussed  above  is  also 
proved  by  the  impact  [energy]  value  listed  in  Table  2 
with  the  material  showing  energy  absorption  in  larger 
quantity  when  compared  to  the  laminated  materials.  It 
also  shows  a  fracture  pattern  similar  to  that  of  the 
chevron  test  specimen,  which  will  henceforth  be  investi¬ 
gated  further. 
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[Text] 

1.  Introduction 

Composite  materials  made  by  laminating  unidirectional 
materials  are  favored  over  ordinary  metals  in  specific 
rigidity  and  specific  strength  and,  in  recent  years,  have 
begun  to  be  used  aggressively  as  components  for  jet 
engines  of  aeroplanes.  Some  practical  applications  in 
static  components  have  been  adready  seen.  In  the  prac¬ 
tical  application  of  FRP  in  a  rotational  component,  e.g., 
a  rotor  blade,  however,  the  damage  produced  to  the 
component  by  colliding  ice,  stones,  birds,  etc.  at  high 
speeds  (Foreign  Object  Damage:  FOD)  poses  a  serious 
problem.  This  report  deals  with  a  basic  experiment  and 
analysis  for  impact  responses  by  use  of  a  cantilever  board 
as  the  model  of  the  blade  in  order  to  assess  the  damages 
produced  in  cases  where  impact  is  applied  to  the  FRP 
blade  at  high  speeds. 

2.  Impact  Experiment 

2. 1  The  Method  of  Experiment 

A  diagram  of  the  test  apparatus  is  presented  in  Figure  1. 

A  steel  ball  (<p  5  mm,  0.5  gm)  was  used  for  the  projectile 
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for  collision  with  the  object  to  be  tested.  The  steel  ball 
was  accelerated  with  gun  powder  to  a  speed  of  100  -  500 
meters  per  second  and  made  to  collide  with  a  stationary 
test  object.  The  speed  of  the  steel  ball  at  collision  was 


Figure  1.  Block  Diagram  of  the  Apparatus 


under  control  by  adjusting  the  amount  of  the  gun  powder 
based  on  some  preliminary  testings.  The  data  obtained 
by  measurements  are  for  speeds  of  the  projectile  at  the 
time  of  collision  and  strain  responses  at  the  time  of 
collision  registered  by  a  strain  ^ge  applied  to  the  test 
object.  The  process  of  deformation  was  recorded  with  a 
high  speed  camera. 

2.2  The  Object  To  Be  Tested 

The  test  object  used  was  a  hybrid  laminated  board 
combining  glass  (G)  and  Kevlar  (K).  The  type  of  lami¬ 
nation  was  plain  weaving  (UD).  TTie  structure  of  the 
lamination  and  results  of  the  three-point  bending  test  for 
stren^h  are  shown  in  Table  1.  The  test  object,  as  a 
simplified  model  for  the  blade,  was  a  flat  board  rectan¬ 
gular  in  shape  and  fastened  in  a  cantilever  state  using 
relevant  tools.  The  shape  of  the  test  piece  and  the  sites  at 
which  the  projectile  collided  and  at  which  the  strain 
gages  were  applied  are  shown  in  Figure  2. 


Table  1.  Laminates 


Type  of  laminate 

A 

B 

C 

Data  for  the  bending  test 

A:  0/90/0 

0°  Rigidity 

90°  Rigidity 

0°  Strength 

90°  Strength 

B:  90/0/90 

G/KUD 

G 

K 

K 

274i 

2446 

87.3 

66.2 

C:  0/90/0/90 

Volume  of 
laminates 
ABABCBABA 

r(so*  ) 


Figure  2.  Simplified  llVlodel  of  the  Blade 


2.3  Results  of  the  Experiments 

The  strain  responses  and  the  damage  mode  produced  at 
the  site  of  collision  by  a  high-speed  collision  test  will 
follow. 

(1)  Damage  Mode 

The  mode  of  the  damage  produced  at  the  time  of 
collision  for  the  FRP  flat  board  is  shown  in  Figure  3, 


VV,  '  '  L 

(a)  Cloth  material  (b)  UD  material 

Figure  3.  Damage  Modd1(opposite  surface) 

along  with  that  for  cloth  material  carried  out  simulta¬ 
neously.  The  speed  of  the  projectile  at  the  time  of 
collision  was  140  meters  per  second,  a  value  below  the 
limit  of  penetration.  As  can  be  seen  from  the  figure,  a 
wide-ranged  separation  of  laminated  layers  was  pro¬ 
duced  in  a  concentric  area  surrounding  the  site  of 
collision  in  the  case  of  cloth  material  and  in  the  direction 
determined  by  the  orientation  of  the  fiber  in  the  case  of 
the  UD  material. 

(2)  Dynamic  Strain  Response 

Figure  4  shows  strain  responses  measured  at  sites  G3  and 
G4  at  the  time  of  collision  with  the  strain  gages  applied 
to  the  test  object  at  these  sites  to  make  measurements. 
The  interval  between  the  adjoining  data  sampling  times 
was  20  p  second  for  Figure  4  (a)  and  0.2  p  second  for 
Figure  4  (b).  It  can  be  seen  from  Figure  4  (a)  that,  after  a 
purturbation  just  after  the  collision  of  the  projectile,  the 
wave  turns  into  a  sine  wave  gradually.  Natural  periods 
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with  frequencies  of  some  1 30  HZ  and  830  HZ  are  also 
observed  in  G3-waves  which  represent  the  frequencies  of 
primary  and  secondary  bendings,  respectively.  In  Figure 
4  (b),  in  turn,  is  displayed  the  wave  on  an  expanded  scale 
for  200  p  second  after  the  collision.  Differences  in  time 
for  the  point  of  the  beginning  of  the  rise  of  the  wave  and 
for  that  of  the  maximum  wave  amplitude  are  notable 
between  sites  G3  and  G4.  For  the  distance  of  the  applied 
gages  of  30  mm,  the  speeds  of  propagation  of  the  waves 
are  computed  to  be  5000  meters  per  second  and  520 
meters  per  second,  and  the  former  value  compares  with 
the  elastic-wave  propagation  speed  where  the  value  for 

{  '^E/P  =  4 
000  i/s) 

Young’s  modulus  in  Table  1  is  used  for  calculation, 
whereas  the  latter  value  may  possibly  represent  the  speed 
of  the  propagation  of  the  wave  for  the  bending  of  the 
board. 


Figure  4.  Strain  Response  to  Impact  (experiment) 


3.  Analysis  of  Impact  Responses 

It  is  possible  to  reproduce  and  assess  the  damage  mode 
providing  the  impact  of  loading  and  the  contact  time  of  the 
projectile  with  the  test  object  at  the  time  of  impact,  among 
other  things,  are  clarified  and  representation  of  the  impact 
response  waves  obtained  in  experiments  by  a  proper  form  in 
analysis  is  possible.  An  analysis  for  impact  responses  using 
the  finite  element  method,  therefore,  was  carried  out,  where 
a  simplified  blade  model  was  used  and  results  compared 
with  those  of  the  above  experiments. 

3.1  The  Method  of  Analysis 

The  general-purpose  FEM-program,  MSC/NASTRAN, 
and  MARC  were  used  for  the  analysis.  Material  param¬ 
eters  in  Table  2  were  used  as  the  shell  elements  suitable 
to  the  analysis  of  laminated  structures. 


Table  2. 

Property 

Lamina  type 

G 

K 

E  1 

4000 

7500 

E2 

900 

600 

G  12 

400 

400 

V  12 

0.26 

0.34 

Before  starting  the  impact  analysis,  analyses  on  the 
natural  frequencies  of  the  blade  model  and  on  relevant 
static  loading  were  carried  out  and  results  obtained  were 
compared  with  relevant  experimental  ones  (Table  3: 
natural  frequency;  Figure  5:  static  analysis)  with  the  view 
to  confirming  the  correctness  of  applying  the  above 
programs  and  the  parameters.  Relevant  analyses  corre¬ 
sponded  well  with  experimental  results  and  sufficient 
precision  was  probable  for  response  analysis.  There  are 
three  types  available  for  the  input  in  the  response  anal¬ 
ysis  as  shown  in  Figure  6:  (a)  The  time  related  shape  of 
the  impact  loading,  (b)  the  initial  velocity,  (c)  the  pro¬ 
jectile  (concentrated  mass  and  its  velocity).  The  (a) 
method  was  the  one  adopted  here. 


Table  3. 


Mode 

KHZ) 

2 

3 

Experi¬ 

ment 

133.5 

345 

835 

Analysis 

MARC 

133.0 

343 

836 

NASTRAN 

136.0..1 

391 

853 

3.2  Results  of  Analysis  for  Impact  Responses 


In  connection  with  low-speed  impacts,  an  analysis  was 
made  in  cases  where  impacts  were  given  to  the  blade 
model  with  an  impact  hammer  incorporating  a  load  cell 
at  its  head.  The  shape  of  the  input  wave,  values  of  the 
loads,  and  contact  time  required  for  the  load  input  had 
been  all  provided  by  relevant  experiments.  Under  this 
condition  the  analysis  was  made  and  the  test  results  were 


JPRS-JST-90-007 
31  January  1990 


13 


Figure  5.  Static-Load  Analysis  (MARC) 
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Figure  6.  Method  of  Input  for  Impact  Analysis 

conflrmed  to  be  same  as  that  of  response  wave-forms. 
(The  results  are  omitted.)  As  for  such  experiments  like 
this  which  involve  high-speed  impact,  however,  it  was 
hard  to  measure  impact  conditions.  Hence  the  value  of 
maximum  impact  load,  one  of  the  input  parameters,  had 
to  be  hypethesized  as  fixed.  Then  by  varying  wave  shape 
and  contact  time,  conditions  best  adapted  to  the  actual 
wave  were  derived.  Changes  of  transient  responses  to 
input  data  are  presented  in  Figure  7.  Final  results  of 
analysis  and  those  of  experiments  for  the  transition  of 
strain  responses  in  the  direction  of  the  span  are  shown  in 
Figure  8,  and  show  good  agreement  as  far  as  the  range 
(within  200  ps)  just  following  the  impact  is  concerned.  In 
particular,  the  time  difference  when  the  strain  wave 
starts  to  rise  and  the  time  difference  of  maximum 
amplitude  between  sites  G3  and  G4  are  4300  meters  per 
second  and  570  meters  per  second  respectively.  These 
figures  with  those  values  were  obtained  by  the  experi¬ 
ments.  The  sine  wave,  as  an  input  wave,  thus  proved 
with  certainty  to  be  capable  of  producing  response  waves 
properly  and  the  contact  time  to  affect  the  simulation 
results  in  a  large  measure.  It  was  also  proved  that  the 
impact  load  affects  the  amplitude  of  the  [response]  wave. 
General-purpose  programs,  meanwhile,  often  offer  a 
procedure  for  impact  analysis  wherein  a  projectile  or 


Wave  Shape 
of  the 


Input  TRANSIEUT  RESPONSE  CG3} 


Figure  7.  Change  in  Response  With  Time-Related  Load 
Patterns  (NASTRAN) 


Figure  8.  Transitional  Process  of  Strain  Response  to  the 
Impact 

concentrated  mass  provided  with  an  initial  velocity  is 
collided  with  the  test  object.  This  method,  which  is 
convenient  enough  to  allow  direct  computation  of  con¬ 
tact  times,  impact  loads,  etc.,  was  used  in  the  present 
analysis  for  the  estimation  of  the  value  of  the  impact 
load  for  reference  sake. 

In  Figure  9  are  shown  the  deformation  modes  in  the 
present  analysis.  It  can  be  seen  that  a  localized  deforma¬ 
tion  propagates  from  the  center,  the  impact  point, 
toward  the  periphery,  with  the  mode  shifting  from  the 
primary  bending  to  the  secondary.  It  was  reported  that, 
for  a  flat  board  fastened  on  all  four  sides,  the  speed  of 
propagation  of  the  maximum  strain-amplitude  agreed 
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Figure  9.  Deformation  Mode  Upon  Impact 

with  the  speed  of  the  progress  of  the  separation  of  the 
laminated  layers  measured  with  the  relevant  high-speed 
photographs. 

3.3  Assessment  of  Damage 

A  final  assessment  of  the  damages  must  be  made  fol¬ 
lowing  the  procedures  so  far  given  above.  Since  com¬ 
posite  materials  show  a  wide  variation  in  strength,  it 
seems  difficult  to  establish  criteria  for  the  quantitative 
assessment  of  the  damages  before  more  data  is  accumu¬ 
lated.  In  the  present  analysis,  therefore,  a  qualitative 
assessment  was  tried  with  inter-layer  shear  stress  and 
Hoffmans’  destruction-rule  formula;  Figure  10  shows  a 
distribution  of  parameters  for  each  of  interlayer  shear 
stress  and  Hoffman  in  a  cross-section  at  the  point  of  the 
impact.  The  parameters  for  interlayer  shear  stresses  grow 
at  the  center  of  lamination  and  those  of  Hoffman  on  the 
back  side  of  the  board  and  particularly  the  fact  that  the 
damage  is  greater  on  the  back  side  rather  than  the  front 
in  the  distribution  of  Hoffman  agrees  with  the  actual 
damage.  It  is  henceforth  necessary  to  verify  various 
destruction  rules  and  to  clarify  the  mechanism  of  inter¬ 
layers  separation  along  with  accumulation  of  data  in 
large  numbers. 

4.  Conclusion 

A  parametrical  analysis  of  experiments  was  conducted 
on  a  model  blade  with  the  view  to  assessing  the  blade 
made  of  FRP,  with  the  following  results. 


Figure  10.  Degrees  of  Damage  at  the  Site  of  Impact 

(1)  An  impact  test  on  a  blade  model  was  conducted  and 
behavior  of  the  blade  assessed  for  the  period  immedi¬ 
ately  following  the  impact  and  for  a  span  that  follows 
that  period. 

(2)  An  impact  analysis  was  made  on  the  basis  of  the 
impact  strain-responses  measured  in  order  to  optimize 
the  load  for  the  impact  and  it  has  been  proved  with 
certainty  that  a  general-purpose  FEM  program  allows 
simulation  analysis  of  composite  materials. 

(3)  The  relevant  parameters  distributed  on  the  basis  of 
Hoffman’s  destruction  rule  proved  that  the  result  agreed 
with  relevant  experimental  results  qualitatively. 
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